Electronic Structure of a Chain-like Compound: TlSe 

§. EUialtioglu,* E. Mete,'^ R. Shaltaf 
Department of Physics, Middle East Technical University, Ankara 06531, Turkey 



o 
o 






I 

o 
o 



> 

oo 

(N 

o 
o 



I 

o 
o 



X 



and 

K. AUakhverdiev,-'^'^ F. Gashimzade,^ M. Nizametdinova,'^ G. Orudzhev^ 

^Institute of Physics, Azerbaijan National Academy of Sciences, Baku 370073, Azerbaijan 

"^Materials and Chemicals Technologies Research Institute, MAM, TUBITAK, Gebze/Kocaeli, Turkey 

^Azerbaijan Technical University, Baku, Azerbaijan 

(Dated: February 2, 2008) 

An ab-initio pseudopotential calculation using density functional theory within the local density 
approximation has been performed to investigate the electronic properties of TlSe which is of chain- 
like crystal geometry. The energy bands and effective masses along high symmetry directions, the 
density of states and valence charge density distributions cut through various planes are presented. 
The results have been discussed in terms of previously existing experimental and theoretical data, 
and comparisons with similar compounds have been made. 

PACS numbers: 71.15.Mb, 71.20.-b, 71.18+y 



I. INTRODUCTION 

Thallium selenide is a III- VI compound with body- 
centered-tetragonal (bet) structure of D^^ {1 4/ mem) 
space groupii The quasi-one-dimensional nature^ of its 
structure makes TlSe a generic example for a series of 
similar binary and ternary chain-like compounds with the 
formula unit of Tl+(TP+Se2~)~ where for the ternary 
compounds monovalent and trivalent cations are of dif- 
ferent atoms like in TlInSe2, TlGaTe2, and TlInTe2. The 
trivalent TP+ ions are surrounded by four tetrahedrally 
bonded Se^^ ions. These tetrahedra share edges to form 
long negatively charged chains of (TP+Sbj")" units that 
are parallel to z-axis coinciding with the optical c-axis 
(see Fig.nj. Monovalent T1+ ions on the other hand are 
surrounded by eight octahedrally positioned Se^~ ions, 
and they electrostatically hold these chains together by 
means of ionic inter-chain forces that are weaker than the 
intra-chain bonds of Tl^+-Se^~ which are ionic-covalent 
in nature. This leads to easy cleavage of TlSe-type crys- 
tals into plane parallel mirror-like plates along c-axis. As 
a result one has a "natural" (110) plane with mirror-like 
surface that is particularly useful for optical measure- 
ments and device applications. Also for the same rea- 
sons these crystals are highly anisotropic in many phys- 
ical properties which consequently are enhanced under 
the influence of high pressures. Mi^ In addition, the spin- 
exchange coupling between T1+ and Tl'^"'" ions is observed 
to be strongeisS. as compared to the intra-chain couplings 
of the same nature between Tl""'' and Tl"+ ions. 

The unit cell contains 8 atoms with atomic positions 
of two T13+ ions at ±(0, 0, c/4), two T1+ ions at ±(a/2, 
0, c/4) and four Se^^ ions at ±(770/2,770/2, c/2) and 
±(770/2,-770/2,0), where 77 is the internal parameter. 
The primitive translation vectors are (— o/2, a/2, c/2), 
(a/2, -a/2, c/2) and (a/2, o/2, -c/2). The lattice con- 
stants for TlSe can be found in the literaturei'^ as a—b— 
8.02±0.01 A and c=7.00±0.02 A with the internal pa- 
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FIG. 1: Crystal structure of TlSe. Thallium atoms are repre- 
sented by larger spheres and selenium atoms by smaller ones 
in order to emphasize the chains for better visualization. Se^~ 
ions are tetrahedrally bonded to TP+ ions forming chains of 
edge-sharing-tetrahedra, along c-axis. The Tl^ ions are the 
spheres with no bonds shown. 



rameter 77= 0.358. 

The covalent bond between trivalent thallium and the 
divalent selenium is of length d(Tl3+-Se2-)= 2.68 A, be- 
ing just a little larger than the sum of the covalent radii, 
1.48 A and 1.16 A, respectively. The length of ionic bond 
between the monovalent thallium placed in an octahe- 
dron of divalent selenium ions is d(Tl"'"-Se^^)= 3.428 A, 
which is smaller than but close to the sum of the respec- 
tive ionic radi:^ of 1.59 A and 1.98 A. Other typical bond 
lengths are d(TP+-Tl"+)= 3.5 A, di{Se'^--Se'^-)= 3.853 
A, d(Tl3+-Tl+)= 4.01 A and d2(Se2--Se2-)= 4.06 A. 

It has been shown that TlSe-type crystals are promis- 
ing materials in device applications as near- and far- 
infrared sensors, pressure sensitive detectors? and as 
7-ray detectors^ Switching phenomenaifi and low- 
temperature metallic conductivity-^- in TlSe have been 



described successfully. 

The energy band gap published by different authors 
are not in accordance and varies between 0.6 to 1.0 eV 
at 300 K}^^ Band structure calculations and compar- 
isons with the existing experimental data showed that 
TlSe-type materials are indirect gap materials; and the 
direct transitions are forbidden according to the symme- 
try selection rulesJi Electronic band structures of TlSe 
and related crystals were obtained by various groups and 
available in the literaturepiii^J^J^iiLiaiia 

Thallium selenide-type crystals posses a three dimen- 
sional electronic nature in spite of their chain-like crystal 
structure. But still the direction normal to the chains 
shows stronger band dispersion and consequently may 
be more conductive. This can be seen from the experi- 
mental values^'^ of direct and indirect gaps for different 
polarization directions. The direct gap is measured as 
0.99 eV when the polarizations of the electric-field vec- 
tor of the incident electromagnetic wave E is normal to 
the optical c-axis (z-axis) whereas when E is parallel to 
c-axis it is measured to be 1.05 eV. Similarly, the indi- 
rect gap values are observed to be 0.68 eV and 0.72 eV, 
respectively for the normal and parallel fields. 

Effect of pressure and temperature on the electronic 
band structure of TlSe was first reported by Gashimzade 
and Orudzhev,'^ and phase transitions in TlSe-type crys- 
tals under pressure were reviewed by AUakhverdiev and 
EUiahioglui^ 

In the present work the results of an ab-initio pseu- 
dopotential calculations using density functional theory 
within the local density approximation for the electronic 
band structure as well as the density of states of TlSe 
are presented. In addition, the valence charge density 
distributions for various atomic planes are calculated. Ef- 
fective masses for various valleys in different symmetry 
directions were estimated by using curvature fit to the 
bands. All these results are compared with the experi- 
mental and other theoretical values available. 



II. METHOD 

We have used a pseudopotential method based on den- 
sity functional theory in the local density approximation. 
The self-consistent norm-conserving pseudopotentials are 
generated by using the TrouUier-Martins scheme^*' which 
is included in the fhi98PP package.^ Plane waves are 
used as a basis set for the electronic wave functions. 
In order to solve the Kohn-Sham equations^^ conjugate 
gradients minimization method^"^ is employed as imple- 
mented by the ABINIT codei^ The exchange-correlation 
effects are considered using the Perdew-Wang scheme^^ 
as parametrized by Ceperley and Alder. ^^ 

Pseudopotentials are generated using the following 
electronic configurations: for Tl, in addition to the true 
valence states (6s and 6p), 5d semicore states are also 
included in the calculation. For Se, 4s semicore and Ap 
true valence states are treated as valence states. The op- 



timized calculation has produced the lattice parameters 
to be a=b= 7.91 A and c= 6.90 A, both of which are 
close to their experimental values within ^^1.4%. 

Good convergence has been obtained for the bulk total 
energy calculation with the choice of a kinetic energy cut- 
off at 20 Ha for TlSe. In the density of states calculations 
the irreducible Brillouin zone (BZ) was sampled with 80 
fc-points using the Monkhorst-Pack^i scheme. 




FIG. 2: The irreducible wedge (heavy lines) of the first Bril- 
louin zone for TlSe structure with the high symmetry points 
and high symmetry lines indicated. 



III. RESULTS AND DISCUSSION 

The first Brillouin zone for TlSe in the bet structure 
is given in Fig. [3 where the high symmetry points and 
high symmetry lines are indicated on the irreducible part 
(l/16th) of the BZ and they are given by: 
R=(0, Tr/a, tt/c), P=(7r/a, n/a, tt/c), N=(7r/a, 7r/a, 0), 
r=(0, 0, 0), T=(0, 27r/a, 0) = T'=(0, 0, 27r/c) and H=(0, 
0, (l-|-u2)7r/c) = H'=(0, 27r/a, (I-u^Jtt/c) where u = c/a. 
Symmetry lines of the Brillouin zone: K={k, 7r/a, tt/c), 
D={Tr/a, n/a, fc), S=ik, k, 0), Z\=(0, fc, 0), A^{0, 0, k) 
and G={k, 27r/a-fc, 0). 

The dashed box in Fig. [21 has edges of \/2{2'K/a) x 
\/2(27r/a) x 2(27r/c), the corners of which are the centers 
(r') of neighboring Brillouin zones. Similarly, the point 
H' is equivalent to H since it is the H-point of the neigh- 
boring BZ. The reciprocal lattice vectors for our choice of 
primitive translation vectors are given by Gi = (0, 27r/a, 
2tt/c), G2=(27r/a, 0, 27r/c) and G3=(27r/a, 2n/a, 0) 
along a*, b* and c* axes, respectively. 

The energy bands calculated for the k-points along the 
high symmetry lines are shown in Fig. |2| Also shown in 
the four rightmost panels are the total density of states 
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FIG. 3: The energy bands for TlSe along the high symmetry lines of the Brillouin zone, the corresponding total density of 
states, and the local densities of states for the Se^~, Tl^+, and T1+, in panels from left to right, respectively. The densities of 
states for the lower valence (semi-core) states due to Tl 5d and Se 4s electrons are shown in a scale ten times smaller than the 
rest of the densities above -9 eV. The top of the valence band is taken to be zero. 



for the TlSe compound, as well as the local densities of 
states for the individual ions, Se^~, Tl^+, and Tl"*", re- 
spectively. 

At the bottom of the figure there are 4 bands originated 
from 4s states of the Se atoms located in the range 13 
to 13.7 eV below the top of the valence band, which is 
chosen as zero. Above this group there are 20 bands 
consisted of 5d states of Tl atoms located in the range 
from -10.6 to -12 eV, and being semicore rf-states they 
are not much dispersed in most fc-directions, except a 
little along D, and along A. 

In the region between -4 and -7 eV, there is an isolated 
group of 4 bands which are made up of mostly the 6s 
states of monovalent Tl ion, and some Se 4p states at 
the bottom of the valence band mixed with the 6s states 
of trivalent cation. Another group consisting 10 bands 
in the upper part of the valence band is mainly due to 
Se 4p-states and 6p states of both Tl atoms. However, 
the uppermost valence band which tops at the symmetry 
point T is composed of mainly nonbonding Se 4p and 6s 
states of monovalent Tl. 



The 2 bands in the lower part of the conduction bands 
are the antibonding mixture of the Se 4p-states with the 
6s states of TP+, and intermixed only slightly (around 
H) with the 12 upper conduction bands that are located 
above ^1.5 eV, and made up of p-states of all three ions. 

All of the bands along the line K joining R and P points 
are doubly degenerate due to time reversal symmetry. — 

The bottom of the conduction band is located almost 
at the midway Di — ^Tr/a, tt/u, tt/2c), along the line D 
joining the points P and N, and corresponds to the ir- 
reducible representation Di. Two additional minima are 
situated along the symmetry line A which connects the 
points r and H, one very close to the midway (A4) and 
the other is very close to the point H. The band at point 
T is a little higher in energy. 

The top of the valence band is sharply defined and lo- 
cated at the high symmetry point T which corresponds 
to the irreducible representation T3. The T3 ^T4 ver- 
tical transition is forbidden in the dipole approximation, 
however, minimal direct transition is allowed at point H 
(near T along A). The bottom of the valence bands is 



also at the symmetry point T with s-like minimum, and 
the valence band width is found to be 7.21 eV. 
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FIG. 4: Energy gap values and effective masses for valleys 
along various directions. The conduction bands are rigidly 
shifted upward in energy by a fixed amount of 0.8 eV (see the 
text) 

The energy gap is underestimated relative to the exper- 
imental value due to the well known artifact of the LDA 
calculationSf2& as a result that the indirect gap appears 
to be negative which instead should be about 0.7 eV. 
Note that the similar deficiency is observed by Okazaki 
et al.~^ for the semiconductor TlGaTe2 using ab-initio 
LAPW method, and resulting to a band structure with 
slightly negative gap also, leading to a semimetal with 
hole pocket at T3 and electron pocket at Di. Hence, for 
a rough correction, if the whole conduction band (of 14 
states) as a block is rigidly shifted upwards by 0.8 eV to 
arrive at the indirect gap value (0.69 eV seen in Fig. ^ 
of TlSe then a direct gap of 1.27 eV is obtained as a 
consequence, which is to be compared with the experi- 
mental value of roughly 1.0 eVf^ The other indirect gaps 
will then be seen as 1.11 eV for the valley along A and 
1.17 eV for the valley along the same direction but just 
a few meV's away from H-point. The curvature of the 
lowest conduction band along TH has the same sign but 
is larger than that of the highest valence band along the 
same edge of the BZ, and thus the vertical gap at H-point 
being 1.22 eV is 0.05 eV smaller than the direct gap at 
T-point where the valence band top is located. 

From the curvatures of bands the effective masses are 
obtained by linear fitting of E versus k^ plots at the close 
proximity of the extrema. They are shown in Fig. ^ for 
different valleys in units of free electron mass TOq- The 
hole effective masses along A and E are found to be 
very close to each other being m*^{T3)=0.33mo, how- 
ever, along TH it is as high as m'^lT^)— 1.16 rrio and in- 
creasing curvature after H along A reduces it to 0.61 TOq 
(at H) giving rise to a combined value of 0.81 mo at T. 
The electron effective masses at T along A and S are 
again similar and have values of to*(T4)=0.37771o and 
0.38 mo, respectively. Along TH the curvature of this 
band is opposite of those along A and G causing T4 to 
be a saddle point as is well-known. Other valleys along 



A correspond to effective masses of 0.34 m-o and 0.22 mo 
as seen in Fig. 0| and the electron effective mass for the 
lowest minimum along D is found to be 0.16 ttIq. The 
experimental values obtained from conductivity and Hall 
effect measurements^^ are TO*=0.3mo and 7TT.^=:0.6m,o 
and from thermoelectric measurements'^*' is Tt,* = 0.86 mo. 




The (001) plane 

FIG. 5: (a) Oblique surface plots and (b) top views of the 
total valence charge density contours for TlSe cut through 
the top (001) plane containing the Se'^" ions (larger balls). 
Tl ions (medium balls) and half of the Se ions (smallest balls) 
are in (004) and (002) planes, respectively. 

The total valence charge densities for different planes of 
atoms were calculated to show the charge transfer which 
are in accordance with the local density of states results 
in identifying the electronic structure of the compound. 
Fig. El shows the plot for the (001) plane which cuts 
through the Se^~ ions at positions ±(770/2,-770/2,0). 
The bottom plane of (002) containing the other Se^^ 
pair of the same tetrahedron gives the same charge den- 
sity plot, except it is flipped in x- or y-axis. 

Fig. El shows the case for the (004) plane that con- 
tains T13+ ions at ±(0, 0, c/4) and T1+ ions at ±(a/2, 
0, c/4). Tl"*" ions having been stripped from their 6p 
electron show their s-like character due to the outermost 
6s^ electrons participating in the valence bands. This 
can also be seen in the energy band picture as rather 
quadratic (s-like) minimum and maximum at the high 
symmetry point T (see Fig. EJ. Donated their 6p and 
6s electrons to bond formation, TP+ ions, on the other 
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The (004) plane 

FIG. 6: (a) Oblique surface plots and (b) top view of the 
total valence charge density contours for TlSe cut through 
the (004) plane containing both Tl^ (center and corners) and 
Tl^+ (edges) ions. Se^~ ions shown are not in plane, but 
either above (larger balls) or below (smaller balls) by c/4. 
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The (220) plane 

FIG. 7: Contour plots of the total valence charge density for 
TlSe cut through the (220) plane containing Tl^+ and Se^" 
ions. Tl^ ions shown are not in plane. The frame shows the 
part of (220) plane inside the cube shown in Fig. |S| 



hand, show some charge density extending towards Se^ 
ions and a negligible amount of d-influence in Fig.|B|[b). 
These findings can only be roughly compared with the 
earlier empirical calculations'^^ since the charge density 
contours for (001), (002) and (004) planes were superim- 
posed in one plot and the amplitude variations along the 
bonds were depicted in a different figure. Therefore, in 
order to compare the pronounced charge accumulation 
at the TP+-Se^~ bonds found in the empirical calcula- 
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FIG. 8: Contour plots of, (a) the total valence charge density 
and (b) the charge density due to the lowest two conduction 
bands, for TlSe cut through an incommensurate plane con- 
taining all three ions. 



tion, charge density distribution on two other planes both 
passing through these bonds are also presented. Fig. [7| 
shows the total valence charge density plots for the (220) 
plane that passes through the TP+-Se^~ bonds, where 
most of the charge is seen to be accumulated on the Se^~ 
ion rather than on the bond. The second plane that 
contains the bonds under consideration contains also the 
monovalent cation. The total valence charge density cal- 
culated for this plane is shown in Fig. Eta), where the 
monovalent cation is seen to be not bonded to the chalco- 
gen ions. Again the TP+-Se^^ bond is seen to be more 
ionic than covalent in nature. Finally, Fig.lHl^b) shows the 
same plane that contains all three ions as (a), however, 
this time the charge density contours depict the contri- 
butions from the lowest two conduction bands only. It is 
clearly seen that the monovalent cation has no electron 
at these antibonding bands which are formed by TP"*" 6s 
and Se^~ Ap states alone, consistent with the rightmost 
panel in Fig. 13 



IV. SUMMARY AND CONCLUSION 

Ab-initio pseudopotential calculations using density 
functional theory within the local density approximation 
has been performed for the first time to investigate the 
electronic properties of TlSe. The energy bands along dif- 
ferent symmetry directions, the local and total densities 
of states, and total valence charge density distributions 
for certain plane-cuts are presented. 

It is shown that the top of the valence band is sharply 
defined and located at the high symmetry point T, which 
corresponds to the irreducible representation T3. The 
bottom of the conduction band is located almost at the 
midway {n/a, 7r/a, 7r/2c) between symmetry points P 
and N, along the the line D, and corresponds to the ir- 
reducible representation Di. 

The distribution, dispersion and the orbital characters 
of the bands are in good agreement with the experimental 
data of the photoemissioni^ The only prior band struc- 
ture calculation for TlSe is reported by Gashimzade et 
al.^^ where they have used a model pseudopotential and 
applied an empirical method. Although there is a general 



agreement between the two results, there are some ma- 
jor differences as well. In their result there is an isolated 
group consisted of the top 2 valence bands which does 
not appear separated from the rest of the valence band 
in our calculations. And the isolated group of 4 bands 
at the bottom of the valence band in our calculations is 
not separated from the rest of the valence band in their 
results. Moreover, the order of their bands at the sym- 
metry point N does not agree for most bands with that 
of the present results. These differences arise from the 
fact that Gashimzade et al>i^ used restricted number of 
plane waves in their empirical pseudopotential method, 
and moreover, the wave functions were not deconvoluted 
enough. It is known that in such approximations it is 
not possible to take into account the screening, and to in- 
clude the electron-ion interaction as well as the exchange- 
correlation effects properly. Another difference with this 
study is, although they have included the 4s semicore 
states of Sc atoms, the 5d states of Tl atoms were not 
taken into account in their calculations. To single out 
the effect of d states, we have made separate calculations 
with and without including the 5d semicore states of Tl. 
Upon inclusion of d states the main change was on the 
4s semicore states of Se which got narrowed down from 
dispersion between -13.7 eV to -12.2 eV to -13.7 eV to 
-13 eV, decreasing the bandwidth from 1.4 eV to 0.7 eV. 
The other change occurred in the lowest two conduction 
bands that are shifted downward and separated from the 
rest of the conduction bands by about 0.2 eV in all sym- 



metry directions. 

It is worth noting here that a similar electronic band 
structure for TlInSe2 is obtained from a calculatiorkiS 
made by constructing the pseudopotentials using the 
scheme suggested by Bachelet et aL'^'^ Although it is a 
semiconductor with a direct gap, valence band structure 
is qualitatively analogous to that of TlSe presented here. 
Isolation of the group of lowest 4 valence bands, the shape 
of the highest valence band, especially along the FTIir 
symmetry directions, and sharpness of the valence band 
top are their similar behavior. Another material to be 
compared is the TlGaTe2, ab-initio band structure of 
which has been published by Okazaki et a\m^ Not men- 
tioning the quantitative differences, in general, the two 
band structures are remarkably similar, apart from some 
variations in the vicinity of T point of the Brillouin zone. 
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